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I.  Introduction  and  Program  Objective 

The  fatigue  initiation  processes  and  subsequent  behavior  of  a  (single)  macrocrack  are  fairly  well 
documented  and  understood  for  conventional,  isotropic  materials.  By  comparison,  the  fatigue 
process  tn  single  crystal  materials  is  significantly  more  complicated. 

Fatigue  damage  in  such  materials  is  observed  to  involve  a  number  of  distinct  'states'.  Further¬ 
more,  transition  from  state  to  state  is  not,  as  it  is  in  more  conventional  materials,  the  simple  pro¬ 
gression  of  accumulated  microdamage  culminating  in  a  fatigue  crack  of  engineering  dimensions. 
In  single  crystals  the  sequence  of  damage  progression  is  greatly  influenced  by  the  regularity  of 
the  microstructue,  its  orientation  with  regard  to  the  principal  loading  axis,  as  well  as  the  more 
obvious  rate-controlling  parameters  such  as  temperature,  AK,  and  state  of  stress.  The  resulting 
macrobehavior  can  therefore  be  somewhat  difficult  to  predict,  as  is  evidenced  by  shifts  in  crack 
propagation  threshold,  multiplicity  of  cracking,  and  changes  in  crack  path  trajectory. 

While  the  fatigue  process  is  commonly  described  as  a  series  of  deterministic  events,  it  is  recog¬ 
nized  that  this  is  influenced  by  a  stochastic  component,  as  is  any  real  physical  process.  This  pro¬ 
babilistic  influence  in  single  crystals  can  be  exploited  to  help  explain  the  macrobehavior. 

This  program  investigates  the  influence  of  microstructure  on  macrobehavior  in  single  crystal  air¬ 
foil  materials  by  considering  the  micromechanics  of  damage  accumulation  as  a  Markov  process. 
The  program  concentrates  on  describing  mathematically  this  relationship  with  the  longer  term 
objective  of  establishing  a  model  to  elucidate  the  behavior  of  this  complex  alloy  system.  It  is 
hoped  that  payoffs  from  this  effort  will  include  suggested  directions  for  alloy  improvement,  as 
well  as  increased  robustness  against  the  influences  of  pernicious  environments  including  hydro¬ 
gen.  Laboratory  data  is  being  obtained  under  separate  funding  so  that  this  program  can  focus  its 
resources  on  the  primary  objective. 

To  summarize,  the  objective  of  this  program  is  to  develop  a  micromechanistic  and  probabilistic 
model  using  the  Markov  paradigm  describing  the  fatigue  process  in  a  class  of  materials  which 
includes  single  crystals  (i.e.,  the  Y  strengthened  superalloys). 


II.  Program  Organization 

The  program  is  structured  into  three  technical  tasks.  The  individual  tasks  are  outlined  here. 
Task  100  -  Micromechanical  Characterization 

This  task  defines  the  mechanisms  of  damage  accumulation  for  the  various  types  of  fracture  ob¬ 
served  in  single  crystal  alloys.  These  fracture  characteristics  will  be  used  to  establish  a  series  of 
Damage  States  which  represent  the  fatigue  damage  process.  The  basis  for  this  investigation  is 
detailed  fractographic  assessment  of  failed  laboratory  specimens  generated  in  concurrent  pro¬ 
grams.  Emphasis  is  on  specifically  identifying  the  micromechanical  damage  mechanisms,  relat¬ 
ing  them  to  a  damage  state,  and  determining  the  conditions  required  to  transition  to  an  alternate 
state.  Principal  investigator  for  this  task  is  Mr.  D.  P.  DeLuca. 
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Task  200  -  Analytical  Parameter  Development 

This  task  will  extend  current  methods  of  fatigue  and  fracture  mechanics  analysis  to  account  for 
microstructural  complexities  inherent  in  single  crystal  alloys.  This  will  be  accomplished  through 
the  development  of  flexible  correlative  parameters  which  can  be  used  to  evaluate  the  crack 
growth  characteristics  of  a  particular  damage  state.  The  analyses  will  consider  the  finite  element 
and  the  hybrid  Surface-Integral  and  Finite  Element  (SAFE)  methods  to  describe  the  microme¬ 
chanics  of  crack  propagation.  Principal  investigators  for  this  task  are  Dr.  S.  E.  Cunningham  and 
Mr.  D.  P.  DeLuca 

Task  300  -  Probabilistic  Modeling 

This  task  will  model  the  accumulation  of  fatigue  damage  in  single  crystal  alloys  as  a  Markov  pro¬ 
cess.  The  probabilities  of  damage  progressing  between  the  damage  states  defined  in  Task  100 
will  be  evaluated  for  input  into  the  Markov  model.  The  relationship  between  these  transition 
probabilities  and  fatigue  life  will  then  be  exploited  to  establish  a  model  with  comprehensive  life 
predictive  capabilities.  Principal  investigators  for  this  task  are  Mr.  C.  G.  Annis,  Jr.  (the  Program 
Manager)  and  Mr.  T.  Watkins,  Jr.. 

III.  Accomplishments 

One  major  accomplishment  is  an  increased  understanding  of  the  micromechanics  of  the  fatigue 
and  fracture  process  in  single  crystal  gas  turbine  blade  superalloys.  We  have  shown  that  the  op¬ 
erative  microscopic  fatigue  crack  fracture  mode  can  affect  the  fatigue  crack  growth  rate.  We 
have  also  shown  that  the  rate  of  input  energy  (a  function  of  stress  intensity  (K),  stress  ratio  (/?), 
and  frequency  (v))  and  test  temperature  (controlling  dislocation  mobility  and  character)  deter¬ 
mine  the  ope  rati,  ■  microscopic  fatigue  crack  fracture  mode. 

To  date  we  have  cataloged  six  distinct  damage  states  based  on  microscopic  fatigue  crack  fracture 
modes.  They  are,  in  order  of  ascending  energy: 

1 .  undamaged  material 

2.  y~Y  decohesion  (submicroscopic  octahedral  fracture  confined  to  the  y  phase) 

3.  microscopic  octahedral  (trans-precipitate)  fracture 

4.  ancillary  decohesion  (may  be  non-crystallographic  fracture  of  the  y  phase)1 

5.  trans-precipitate  non-crystallographic  (monoplanar)  fracture 

6.  trans-precipitate  non-crystallographic  (ancillary)  fracture 


These  discrete  states  of  fracture  energy  dissipation  (separated  by  transition  regions)  are  controlled 
by  dislocation  dynamics.  They  can  possibly  be  modeled  by  Arrhenius-like  rate  kinetics  to  obtain 
an  apparent  activation  energy  for  each  discrete  transition  event  (or  end  point,  i.e.,  and  *IC). 


1  This  mode  has  been  paitially  documented,  but,  as  yet  is  not  fully  understood.  A  temperature  gradient  test 

between  room  temperature  and  800°F  is  planned. 
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We  hypothesize  that  these  activation  energies  then  set  the  energy  limits  between  which  a  particu¬ 
lar  fracture  mode  (with  its  characteristic  da/dN  behavior)  will  be  operative. 

Markov  Paradigm 

•  From  the  perspective  of  the  Markov  Paradigm  these  micromechanical  events  represent  the 
transitional  and  absorbing  states  that  determine  the  physical  outcome  of  the  fatigue  process.  The 
Markov  model  is  described  in  greater  detail  in  the  attached  ASME  paper,  "Markov  Fatigue  in 
Single  Crystal  Airfoils"  by  Charles  Annis  and  Daniel  P.  DeLuca,  presented  at  the  International 
Gas  Turbine  and  Aeroengine  Congress  and  Exposition,  Cologne,  Germany  June  1-4,  1992. 

•  We  have  now  identified  what  appears  to  be  the  bulk  of  the  microscopic  fracture  modes. 
Significance 

The  identification  and  cataloging  of  these  fracture  states  along  with  the  correlation  of  fatigue 
crack  growth  (FCG)  behavior  represents  a  significant  contribution  to  the  study  of  fatigue  and 
fracture  and  thus  increased  fatigue  resistance  through  microstructural  modifications  based  on 
those  results.  For  more  details  on  microstructural  modifications  which  have  resulted  in  an  order 
of  magnitude  increase  in  LCF  life  in  hydrogen  and  a  50%  reduction  in  fatigue  crack  growth  rate 
see  our  contract  progress  reports  FR-21998-02  (15  December,  1991)  and  FR-21998-08  (15  June, 
1992);  also  see  "Improved  Crack  Growth  in  Hydrogen  With  Modified  Precipitate  Morphology 
Single  Crystal  Nickel",  by  D.  P.  DeLuca  et  al..  Proceedings  of  the  Conference  on  Advanced 
Earth-to-Orbit  Propulsion  Technology  1992. 

Planned  Efforts  Remaining 

We  will  quantify  Markov  transition  probabilities  based  both  on  statistical  techniques,  such  as 
maximum  likelihood  parameter  estimation,  and  through  quantitative  metallography  identifying 
the  relative  proportion  of  damage  in  the  neighborhood  of  the  advancing  crack. 

We  will  continue  to  correlate  fracture  details  with  experimental  FCG  data  currently  being  devel¬ 
oped  under  a  related  program.  We  will  also  be  analyzing  temperature  gradient  FCG  tests  and 
fractures  to  focus  on  transition  behavior  for  better  resolution  of  these  events.  We  will  commence 
Task  200,  Analytical  Parameter  Development,  and  continue  Task  300,  Probabilistic  modeling. 

Presentations  Given  and  Papers  Published 

The  paper  "Markov  Fatigue  in  Single  Crystal  Airfoils,"  by  Charles  Annis  and  Daniel  P.  DeLuca 
was  presented  at  the  International  Gas  Turbine  and  Aeroengine  Congress  and  Exposition,  Co¬ 
logne,  Germany  June  1-4,  1992. 

The  paper  "Improved  Crack  Growth  in  Hydrogen  With  Modified  Precipitate  Morphology  Single 
Crystal  Nickel",  by  D.  P.  DeLuca,  H.  B.  Jones,  B.  A.  Cowles,  and  F.  D.  Cobia  was  presented  at 
the  Conference  on  Advanced  Earth-to-Orbit  Propulsion  Technology  held  at  NASA  George  C. 
Marshall  Space  Flight  Center,  Huntsville,  Alabama,  1992. 
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This  program  investigates  the  seemingly  unusual  behavior  of  single  crystal  airfoil 
materials.  The  fatigue  initiation  processes  in  single  crystal  (SC)  materials  are 
significantly  more  complicated  and  involved  than  fatigue  initiation  and  subsequent 
behavior  of  a  (single)  macrocrack  in  conventional,  isotropic,  materials.  To 
understand  these  differences  is  the  major  goal  of  this  project. 
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Markov  Fatigue  in  Single  Crystal  Airfoils 

CHARLES  ANNIS  and  DANIEL  P.  De  LUCA 

Pratt  &  Whitney 

West  Palm  Beach,  Florida  33410-9600 


Abstract: 

This  paper  considers  the  influence  of  microstructurc  on 
macrobehavior  in  single  crystal  airfoils  bv  treating  the 
micromcchnnics  of  damage  accumulation  as  a  Markov 
process.  Single  Crystal  Fatigue  (SCF)  is  a  result  of  several, 
simultaneous  (competing),  damage  mechanisms,  which  arc 
selectively  favored  by  particular  combinations  of  external 
conditions.  As  with  any  real  physical  process.  SCF  is  also 
influenced  by  a  stochastic  component.  This  probabilistic 
influence  can  be  exploited  to  help  explain  the  macrobchavior. 

We  begin  with  a  description  of  single  crystal  materials  and 
how  they  differ  from  more  conventional  (isotropic)  alloys. 
Relationships  arc  suggested  among  the  more  probable  of 
several  competing  microstructural  damage  mechanisms  and 
specific  rate-controlling  parameters.  The  states  of 
microstructural  damage  are  then  described  and  catalogued, 
and  the  various  avenues  of  damage  accumulation  arc 
investigated.  Next,  the  Markov  paradigm  is  reviewed  as  it 
applies  to  these  materials.  Finally,  a  Markov  model  is 
presented  to  describe  the  rather  complex  behavior  observed  in 
single  crystals,  and  its  use  in  liflng  gas  turbine  engine  airfoil* 
is  discussed. 


Introduction 

The  difficulty  of  fatigue  life  prediction  for  jet  engine  airfoils 
is  cxaccrbaled  by  the  complex  behavior  of  single  crystal 
materials.  While  the  fatigue  initiation  processes  and 
subsequent  behavior  of  a  (single)  macrocrack  ate  fairly  well 
understood  and  documented  for  conventional  isotropic 
materials,  fatigue  in  single  crystal  metals  is  signifcantlv  more 
intricate. 

Fatigue  in  such  materials  is  observed  to  involve  a  small 
number  of  distinct  states.  These  damage  states  are  mt_  simply 
numerous,  successively  increasing  lengths  of  i  single 
macrocrack,  but  rathe;  microstructiirally  discernible 
conditions  which  comprise  the  damage  accumulation  process 


itself.  Furthermore,  transition  from  state  to  slate  is  not.  as  it 
is  in  more  conventional  materials,  the  repeated  progression  of 
accumulated  microdamagc  culminating  in  an  incremental 
advance  of  the  fatigue  crack  front.  In  single  crystals  t he 
sequence  of  damage  progression  is  greatly  influenced  by  the 
regularity  ot  the  microstructurc.  its  orientation  with  regard  to 
the  principal  loading  axis,  as  well  as  the  more  obvious 
rate-controlling  parameters  such  as  temperature,  loading 
frequency,  and  state  of  stress.  The  resulting  m.  cn  behavior 
can  therefore  be  somewhat  unexpected,  as  is  c>  idenced  by- 
shifts  in  crack  propagation  threshold  (Telcsman  and  flhosn 
1 989).  multiplicity  of  cracking,  and  unexpected  changes  in 
crack  path  trajectory  (DcLnca.  ct.nl.,  1991V 

C  haracteristics  of  Single  Crystal  Materials 

Modern  gas  turbine  flight  propulsion  systems  employ  single 
crystal  materials  for  turbine  airfoil  applications  because  of 
their  superior  performance  in  resisting  creep,  oxidation,  and 
thermal  mechanical  fatigue  (TMF).  These  properties  have 
been  achieved  by  composition  and  alloying,  of  course,  but  also 
by  appropriate  crystal  orientation  and  associated  anisotropy. 

F.arly  aeroengine  turbine  blade  and  vane  materials  were 
conventionally  cast,  equiaxed  alloys,  such  as  INI00  and 
Rene  KO.  I  his  changed  in  the  Into  I9MK  with  the  irmoduction 
of  directionally-solidified  (DS)  MAR-M200  <  Hf  airfoils,  flic 
DS  process  produces  a  <001  >  crystallographic  orientation, 
which  in  supcralloys  exhibits  excellent  strain  controlled 
fatigue  resistance  due  to  its  low  clastic  modulus.  T  in  absence 
of  transverse  grain  boundaries,  r  60%  redumion  in 
longitudinal  modulus  compared  with  equiaxed  grains,  and  its 
corrcsponding  improved  resistance  to  thermal  fatigue  and 
creep,  permitted  significant  increases  in  allowable  metal 
tcmpcratuics  and  blade  stresses.  St  ill  further  progress  ws* 
achieved  in  the  mid  1 970s  with  the  development  of  single 
crystal  airfoils  (Gcll.  et.al..  1980). 

The  First  such  material,  PWA  MHO,  has  a  considerably 
simpler  composition  than  preceding  cast  nickel  blade  alloys 
because,  in  the  absence  of  grain  boundaries,  no  grain 


Presented  et  the  Intemetlonal  Gas  Turbine  end  Aeroengine  Congress  end  Exposition 
Cologne.  Germany  June  1-4,  1992 


boundary  strengthening  elements  arc  required.  Deleting  these 
grain  boundary  strengthened,  which  are  also  melting  point 
depressants,  increased  the  incipient  melt  temperature.  This, 
in  turn,  allowed  ncatiy  complete  y'  solutioning  during  heat 
treatment  and  thus  a  reduction  tn  dendritic  segregation  The 
absence  of  grain  boundaries,  the  opportunity  for  full  solution 
heat  treatment,  and  the  minimal  post-heat  treat  dendritic 
segregation,  result  in  significantly  improved  properties  as 
compared  with  conventionally  cast  ot  directionally  solidified 
alloys.  Single  crystal  castings  also  share  with  Its  alloys  the 
'001  ''  crystal  orientation,  along  with  the  benefits  of  the 
resulting  low  modulus  in  the  longitudinal  direction. 

brail  Whitnev  has  developed  numerous  sinelc  crvstal 
materials.  Like  most,  P\Y  A  I4S0  and  PW  A  1 4 S4  ate  , 
strengthened  cast  mono-grain  nickel  supers  Hoys  based  on  the 
Ni-Cr-AI  system.  The  bulk  of  the  microstructurc  consists  of 
approximately  r»f)n <>  by  volume  of  cuboidal  ,  ’  precipitates  in 
a  y  matrix.  The  precipitate  ranges  from  0.35  to  0.5  microns 
and  is  an  ordered  Face  Centered  Cubic  (FCC)  nickel 
aluminidc  compound.  The  macrostructurc  of  these  materials 
is  characterized  by  parallel  continuous  primary  dendrites 
spanning  the  easting  without  interruption  in  the  direction  of 
solidification.  Secondary  dendrite  arms  (perpendicular  to 
solidification)  define  the  inicrdcndritic  spacing.  Solidification 
lor  both  primary  and  secondary  dendrite  arms  piocccds  in 
()0I  '  type  crystallographic  directions.  I  ^dissolved 
eutectic  pools  and  associated  microporosity  reside  throughout 
the  inicrdcndritic  areas.  These  features  act  as  microstructura' 
discontinuities,  and  often  exert  a  controlling  influence  on  the 
fatigue  initiation  behavior  of  the  alloy.  Also,  since  the 
eutectics  arc  structurally  dissimilar  from  the  surrounding 
matrix  their  fracture  characteristics  will  differ. 

Single  Crystal  Fatigue 

The  fatigue  process  in  single  crystal  airfoil  materials  is  a 
remarkably  complex  and  interesting  process.  In  east  single 
crystal  nickel  allovs.  two  basic  fracture  modes, 
crystallographic  and  non-erystallographic.  are  seen  in 
combination.  The'  occur  in  varying  pioportions  depending 
upon  temperature  anil  stress  state.  ( Tvstallographic 
orientation  with  respect  to  applied  load  also  a  fleets  the 
proportion  of  each  and  influences  the  specific  crvstnliogrnphic 
planes  and  slip  dilations  involved.  Mixed  mode  fracture  is 
observed  under  monotonic  as  well  as  cy  clic  londi’mrs. 

'single  crystal  turbine  blades  ate  cast  such  that  the  radial  axis 
of  the  component  is  essentially  coincident  with  the  001 
crystallographic  direction  which  is  the  direction  of 

solidification.  Crystallographic  fracture  is  usually  seen  a1 
either  octahcdial  alone  multiple  (III)  planes  or  under  certain 
circumstances  as  (00|i  cleavage  alone  cubic  planes  I  he 
fatigue  and  tincture  processes  can  be  visualized  by  relating 
the  frai  ioeraphic  lea  tines  to  structtnal  aspects  of  the  allov 
system.  In  Figure  I  the  FCC  unit  teT  structure  is  shown 
along  with  a  depiction  of  the  (III)  crystallographic  plane 
Such  planes  in  eight  neighbor ine  unit  cclis  describe  an 
oetahedion  The  tractate  ot  an  I  Cl  specimen  shown  in 
Figure  2  displavs  ptomment  (III)  octahedral  s'ne  n  plane 

Non-cry  sialltrgtaphie  traeture  is  also  observed.  low 
temperatures  lain:  et  v  stallographu  li.utute  V  htghet 
temperatures,  m  i  m  ranee.  ■,n;;ll  amounts  ot 

non-crystalloeraphii  ptopagatton  have  the  appearance  ot 
t  ran  sura  tn  i  Ft  tattuue  in  a  related  One  twain  ouuinv.v  allov 
Lnder  some  londutons.  ' ins  propagation  ehantv.  aireos' 
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immediately  to  the  highly  crystallographic  mode  along  till) 
shear  planes,  frequently  exhibitinij  prominent  smations 
emanating  from  the  fatigue  origin  and  continuing  to  the  final 
fracture  legion.  Under,  other  conditions  the  nun- 
crystallographic  behavior  can  continue  until  global  failure 
occurs,  as  is  seen  in  the  turbine  blade  fracurc  surface  shown 
in  Figure  V  At  intermediate  temperatures  (around  TMC] 
non*':ryst allograph!*'  propagation  is  more  pronounced  and 
mav  continue  until  tensile  overload,  or  may  transition  tr 
stihcritic.il  crystallographic  propagation  as  at  427C.  At  9X2C 
propagation  is  almost  entirely  non-crystallouraphic.  similar  to 
transeranul.u  propagation  in  a  polycrystal. 


rigurc  V  Nnn-enstallngraphic  Irnnsprccipitalc  fracture  of  a  turhinc 
Mailt1:  Siii’li  failure  t  an  occur  normal  to  lilt'  print  ipai  loading 
,»\is  al  rlcsaied  temperatures  ami  lower  stresses. 


Damage  (  atalogue 

Del  uca.  ei.al.  ||99|)  have  compiled  descriptions  of  single 
crystal  fracture  morphologies  and  have  postulated  fatigue 
mechanisms  for  them.  (  This  work  is  continuing  and  no  claim 
is  made  that  these  are  the  only  mechanisms  by  which  single 
crvstnl  ma'enals  can  experience  failure.)  It  is  helpful  to 
consider  these  different  mechanisms  as  unique  damage  states 
They  air 

I.  undamaged  material 


3.  dislocation  penetration  of  /  on  limited  (III)  planes 

4.  local  octahedral  free  surface  formation 

A  dislocation  penetration  of  on  additional  slip  systems 

<>.  cross  slip,  dislocation  climb,  and  precipitate  bypass:  creep 

7.  global  crystallographic  octahedral  fracture 

X.  non-n ystallographic  iransprccipitate  fracture 

The  first  six  of  these  arc  'ransitum  whites,  the  final  two  arc 
absorbing  states.  The  distinction  is  discussed  in  later 
paragraphs. 


Considerable  progress  has  been  made  in  understanding 
micromechanical  causc-and-cfTect  relationships  (Drcshficld. 
1986;  Tclcsman  and  Ghosn.  1989;  DeLuca  ct.al..  1991).  With 
increased  understanding,  however,  comes  the  realization  that 
many  of  the  microstructural  causes  can  only  be  determined 
with  electron  microscopy  and  other  e  x  post  facto  methods. 
Therefore,  a  fatigue  life  prediction  system  requires,  in 
addition  to  an  understanding  of  the  physical  phenomena,  a 
method  for  cataloging  the  relative  likelihoods  of  these 
micro-occurrences,  and  a  formalism  for  describing  their  likely 
interactions,  the  collective  behavior  of  which  constitutes  ihe 
observable  fatigue  process.  The  Markov  paradigm  seems 

ideally  suited  to  this  task.  Put  another  way.  micromcchanics 
describes  the  physical  relationships  linking  damage 
accumulation  and  microstructurc.  But  because  fine 
microstructurc  is  random,  this  information  cannot  be  used 
a  priori.  Thus,  a  stochastic  model  is  useful  for  damage 
prediction. 


Markov  Fatigue 

The  fatigue  and  fiacturc  processes  in  single  crystal  materials 
can  be  viewed  as  a  sequence  of  micromechanical  events. 
Some  of  these  events  are  transitional  in  that  they  can  bo 
emerged  from  to  enter  a  more  damaged  state:  therefore, 
damage  is  accumulated  in  the  transition  from  one  state  to 
another.  Other  states  arc  non-transitional.  or  absorbing 
states,  representing  ultimate  material  degradation  such  as  the 
creation  of  a  thermodynamic  free  surface,  or  crack.  (Since 
fatigue  is  essentially  irreversible,  transition  probabilities  from 
a  more  damaged  state  to  one  of  lesser  damage  are  defined  ns 
zero.) 

Because  ( besc  processes  .it  furulamcnttilly  probabilistic  in 
nature,  the  entire  ensemble  can  be  modeled  as  a  Markov 
process.  Here  the  aggregate  (macroscopic)  material  response 
can  be  described  as  having  resulted  from  sequences  of 
microevcnts.  While  the  underlying  model  enjoys  a  certain 
simplicity,  it  is  remarkably  powerful  in  describing  complex 
interactive  behavior. 

Consider  an  elemental  volume  of  material  being  subjected  to 
cyclic  loading  at  some  operating  temperature,  for  example  a 
fatigue  specimen  gage  volume.  At  any  given  cycle  count,  n. 
(or  time,  or  block  of  cycles)  the  damage  process  is  in  one  of 
the  states.  1.  2.  3. ...  A.  Let  ,Y„  represent  the  damage  process 
at  time  n.  so  that  the  observed  sequence  of  states  is 

.V,,  ,Y; . Y„,  Y„ .  and  so  on.  Under  certain  conditions.’ 

this  sequence,  or  chain  of  cvenis  fa  Markov  chain)  can  be 
summarized  by  an  A  x  A  matrix  of  transition  probabilities.  P 
(equation  I). 

Now  let  />..  be  the  probability  of  transition  from  state  /  to  state 
j.  Hus  is  the  conditional  probability  trust  the  fatigue  process 
will  be  in  state  ,Y„  =  j  at  time  n.  given  I  hat  the  preceding  state 
was  ,Y„  ;  —  i.  that  is  p,,  -  Prnh[X„  j\  ,Y„  ,  -  i  I.  If  the 
damage  process  is  at  state  i  at  time  n  I.  then  il  must  be  in 
some  other  (but  not-ncccssarily  different)  state  at  time  n.  The 


1  Fhc  Markov  here  are  c^cnhaJIv  that 

(1)  1  lie  future  depends  only  upon  the  present,  and  iioi  upon  the  past  Staled 
differently:  the  probabilities  of  future  events  depend  only  on  the  current  stale 
the  path  which  load  to  that  stale  tc  irrelevant.  Symbolically  fhi«  k 

Proh[  V,  -  i  |  Vi  -  / 1  .  Vj  =*  i}  .  ...  Y„  i  - 
-  Pnl rf.V.  -/|  V.  ,  -  / 

(2)  Ihcsc  probabilities  arc  constant  over  time  I  hi*  provide*  the  theoretical 
justification  for  equation  V 
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movement  from  state  i  to  state  /  is  controlled  be  the  matrix 
of  conditional  probabilities.  I 'sing  X  damage  states  ns  an 
example: 


r  i  p  ■  ■■  r  -• 

/>::  ■■■  />:« 

P  =  Eqn  I 

/>■,  />';  ...  p-S 


I  ite  state  /  is  called  an  absorbing  state  if  it  cannot  be  departed 
once  it  has  been  entered.  Thus,  if  i  is  an  absorbing  state,  then 
p  -  I  and  p.  -  0  for  all  i  and  /  -r  i.  It  is  important  to  notice 
'lint  there  mnv  be  more  than  one  absorbing  state,  i . e . :  more 
'han  one  type  of  failure,  for  example,  global  crystallographic 
octahedral  fracture,  or  failure  b'  non-cryxi, biographic 
transprccipitatc  fracture. 

This  statistical  groundwork  permits  answers  to  two  main 
questions: 

■  What  is  the  expected  time,  or  number  of  accumulated 
cycles,  required  to  produce  failure  from  a  given  state'.’ 

■  What  is  the  probability  of  an  early  failure? 

These  questions  can  be  addressed  in  terms  of  the  transition 
probabilities.  Here  is  a  greatly  simplified  example. 


Transition  Probabilities  and  Fatigue  Lifetimes 

There  is  a  relationship  between  transition  probabilities  and 
fatigue  lifetime.  Consider  the  simplest  situation  with  only  two 
slates:  state  I  is  damaged,  but  unfracturcd.  and  state  2  is 
fracture.  Obviously,  state  2  is  an  absorbing  state.  Sow  ft,,  is. 
by  definition,  the  probability  of  returning  to  state  I.  at  cycle 
eount  n.  from  already  being  in  state  1  at  cycle  count  n  -  I. 
Similarly,  p,:  is  the  probability  of  becoming  fractured  (state 
2)  at  cycle  count  ».  given  being  unfracturcd  (state  I)  at  cycle 
count  ii  -  I .  Clearly,  p2,  —  0  and  pn  -  I.  that  is.  the 
probabilit v  of  becoming  whole  after  experiencing  fracture  is 
zero,  and  the  probability  of  remaining  fractured  is  one 
Notice,  too,  that  />,,  4-  p,2  -  I.  since  if  the  system  is 
currently  in  stale  I  Mhc  first  subscript)  then  it  must  be  in 
either  state  1  oi  si  ne  2  (the  second  subscript)  at  the  next 
interval;  there  arc  no  other  alternatives.  Thus  the  sum  of  the 
probabilities  must  be  100%.  Now,  if  /»,,  is  large  (i.e.  nearly 
one),  then  />, 2  is  necessarily  small  (ncatlv  zero),  and  the 
associated  cyclic  life  should  he  longer  than  for  a  situation 
where,  sa  v.  p,,  -  p  ,  -  1/2. 

Tor  this  gross  over-simplification,  consider  a  sequence  of 
events  culminating  in  failure: 

cycle:  I.  2.  T  4 .  n  I .  n 

state:  I.  I.  I,  I .  I.  2 


lilt)  —  p,i  *  p,,  X  p.;  x  p„  x  ...  x  />,,  y  p,~  =  p.'ll  -  p  -)" 
.  r»'pra*c«1  n  -  ' 

F.qn  2 

which  is  the  probability  density  of  a  Geometric  distribution, 
sometimes  referred  to  as  the  waiting  time  distribution. 

For  this  distribution  the  expected,  or  average .  ;•  lies  to  failure 
can  be  shown  to  be  n  —  \lp,2  (Mood.  G’-iybilt.  anil  Rocs. 
|974).  For  example,  if  p„  =  0.99  then  pti  ~  0.01  and 
n  —  100.  That  is.  if  the  probability  of  failure  an  any  given 
evclc  is  one  in  one  hundred,  then  the  average  cyclic  life  is  100 
cycles  to  failure.  If  p„  =  0.999,  then  p,2  —  0.001  and 
il  —  1/0. 001  —  1000  cycles.  The  point  is  that  cvclic  life  is  a 
natural  outcome  of  the  state  transition  probabilities. 

Notice,  too.  that  the  properties  of  the  geometric  distribution 
can  exploited  to  estimate  the  probability  of  an  early  failure, 
where  n  •  '  /7.  In  this  very  simple  example,  the  probability 
of  a  failure  before  A  counts  (or  cycles)  is 

P{n'  V)  =  T  pi  I  -  pV  1  Eqn  3 

r  I 

where  p  is  the  probability  of  failure  on  the  next  cycle,  p in 
this  example. 

The  situation  is  analogous  for  the  more  realistic  model  in 
which  fatigue  damage  is  catalogued  into  more  than  two 
damage  states.  The  arithmetic  and  subsequent  statistics  arc 
considerably  more  involved,  however. 

It  can  be  shown  (Taylor  and  Karlin.  1984)  that  by 
partitioning  the  transition  matrix.  Equation  1.  into 
transitional  and  absorbing  portions. 


P- 


Q  I 

— i — 


n  I  / 


Eqn  4 


with  r  transitional  stnics  anil  <;  absorbing  states,  then  the 
t  y  t  matrix  O  is  the  portion  of  P  corresponding  to  transitional 
(non-absorbing)  states,  and  V  is  a  t  -  a  matrix  representing 
the  probabilities  of  moving  from  a  transitional  state  to  an 
absorbing  state.  /  is  an  identity  matrix  of  appropriate 
dimensions. 

It  can  be  further  shown  that 

■  /.  =  (/  -<?)'!  Eqn  5 

where  I  is  a  column  matrix  of  ones,  a, id  L  is  the  vector 
representing  times  to  failure,  from  a  given  state. 

■  The  expected  time  to  failure  of  the  undamaged  material 
is  the  first  clement  of  this  matrix.  L,. 

Thus  the  desired  engineering  pninmctcr.  lifetime,  is  obtained 
from  the  Mai  kov  description  of  t he  fatigue  process. 


The  probabilit  v  of  (exactly)  n  cycles  to  failure  is  then. 
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Istimating  the  Model  Parameters 

Consider  again  the  elements  comprising  the  transition 
probability  matrix  in  Initiation  I.  It  is  necessary  that  these 
elements  he  estimable  tor  the  Markov  model  to  he  a  useful 
engineering  tool.  Some  elements  arc  easily  determined. 
Because  fatigue  is  irreversible,  transitions  from  a  more  to  a 
lesser  damaged  slate  nccessarilv  have  probability  zero.  Thus, 
ill  the  probabilities  heiow  the  matrix  diagonal  in  liquation  I 
are  zero,  ie:  />.  -  0,  for  /  /  For  this  material  there  appear  to 

he  two  absorbing  states,  global  crystallographic  octahedral 
fracture  (state  and  failure  by  non-crystallographic 
transprecipitate  fracture  (state  8).  Tints  n-~  I  and  /u.  I 
so  that  the  lower  right  portion  of  the  matrix  is  the  2  ■  2 
identity  matrix.  Now.  since  />••  -  1  and  /»«  —  I.  and  because 
the  sum  of  the  probabilities  within  a  row  must  be  unity,  all 
other  entries  in  rows  7  and  8  must  be  zero.  Since  there  arc  six 
transition  states  and  two  absorbing  states,  the  partitioned 
matrix  is  then. 


D,  . 
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The  remaining  transition  probabilities  van  be  estimated  from 
experimental  observations,  which  include  the  relative 
proportions  of  crystallographic  octahedral  fi.ieture  and 
non-crystallographic  transprccipitatc  fracture,  and  their 
average  cyclic  lives.  The  computations  arc  conceptually 
straightforward  but  tedious  in  practice.  The  idea  is  to  invoke 
the  familiar  least-squares  criterion  and  select  ns  probability 
estimates  those  values  which  minimize  the  summed  squared 
errot  between  these  observed  quantities  and  those  predicted 
bv  the  model. 


Conclusions 

Pratt  Whitncv  has  enjoyed  modest  success  in  using  this 
Markov  model  to  explain  some  of  the  unusual  behavior 
observed  in  single  crystal  airfoil  materials  in  fatigue. 
Although  only  in  its  early  Pcvclopmcni.  ii  promises  to  provide 
a  workable  life  prediction  system  for  these  complex  materials, 
something  which  has  until  now  proved  rather  elusive.  The 
method  is  also  being  evaluated  for  describing  damage 
accumulation  in  composite  materials,  again  where  damage  b 
complex  and  docs  not  lend  itself  to  description  in  terms  of  a 
single  physical  entity,  such  as  cracklcngth. 
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